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We have performed cascade genetic algorithm and ab initio atomistic thermodynamics under the 
framework of first-principles density functional theory to study the (meta-)stability of a wide range 
of LiajNj, clusters. We found that hybrid xc-functional is essential to address this problem as a 
local/semi-local functional simply fails even to predict a qualitative prediction. Most importantly, 
we find that though in bulk Lithium Nitride, Li rich phase, i.e. LiaN, is the stable stoichiometry, 
in small Lia,Ny clusters N-rich phases are more stable at thermodynamic equilibrium. We further 
show a that these N-rich clusters are promising hydrogen storage material because of their easy 
adsorption and desorption ability at respectively low (< 300K) and moderately high temperature 
(> 600K). 
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Hydrogen is the most abundant element in the universe 
that contains the highest energy density per unit mass, 
and when it burns it produces only water and energy. 
In view of this, it has been considered to be one of the 
promising solutions for the clean alternative energies Il¬ 
ls]. Significant amount of research initiatives have been 
taken place along this direction mz], but unfortunately, 
hydrogen is needed to be produced from water or other 
solid-state materials as it is not freely available in nature. 
Therefore, although hydrogen can be considered as re¬ 
newable, if it is produced from water, it costs more energy 
to get produced, than the energy one recovers on burning 
it. This necessitates, designing of promising solid-state 
materials or nano-structures, where hydrogen molecules 
can be stored. This has been, therefore, an active field 
of research for the last few decades and various kind of 
materials such as metal hydrides IIIHIE], clusters [la¬ 
in], nanostructures [THH21], highly porous metal-organic 
frameworks (MOFs) [51151 E5IB5] . clathrate hydrates [251 - 
E], covalent organic frameworks (COFs) |32J [33] etc. 
have been studied extensively. In particular, complex 
binary hydrides involving light metals such as Li, Mg, 
Ca etc. have been extensively investigated because 
of their high gravimetric storage capacity. Amongst 
them, lithium nitride [LiaN], lithium imide [Li 2 NH] or 
amide[LiNH 2 ] have been found to exhibit strong affinity 
for H 2 because of the reversible reactions Li 3 N-|- 2 H 2 o 
Li 2 NH-bLiH-bH 2 O LiNH 2 -b2LiH [34]. Unfortunately, 
though this material can store sufficiently high amount 
of hydrogen, the desorption kinetics is poor for its on¬ 
board practical application. Li atoms being ionized as 
Li-*- cations, strongly attracts [NH 2 ]“ complex anion and 
thereby form a strong bond, making the desorption of 
hydrogen to be high |7]. In order to lower the desorp¬ 
tion temperature to a permissible range, doping suitable 
metal ad-atoms have already been tried [T] |3S] . 

We have explored the possibility of realising clusters of 
these materials. It has been found that reducing the 


number of particles in a cluster reveals the possibility 
of several interesting properties. In a range where mat¬ 
ter is reduced to sizes of only a few atoms, the intrin¬ 
sic properties of the so-called clusters are non-scalable 
from their bulk analogues. We expect that significant 
amount of hydrogen adsorption would be possible in a 
non-dissociative manner due to presence of high surface- 
to-bulk ratio and due to the weak dispersive interactions 
the hydrogen desorption could be favourable. In past, 
it’s been studied to understand the hydrogen storage ef¬ 
ficiency of small lithium amide clusters [LiNH 2 ]n using ab 
initio molecular orbital theory |36j . But materials prop¬ 
erty changes under operational environment (e.g. tem¬ 
perature (T) and pressure (p) in an atmosphere of re¬ 
active molecules). The thermodynamics and kinetics at 
the relevant temperature (T) and the nature of the en¬ 
vironment determine the composition and structure of 
clusters. In thermodynamic equilibrium, only structures 
and compositions that minimize the Gibbs’ free energy of 
formation of the composite cluster-l-ligands(gas) system 
will be the most stable. Therefore, one has to first ensure 
the most stable phases and compositions of such clusters 
at thermodynamic equilibrium. And after that the ef¬ 
ficiency of hydrogen storage should be aimed at. Till 
date, no such consistent theoretical study has been per¬ 
formed, in the field of hydrogen storage in small clusters, 
that has considered the (meta-)stability of the clusters 
by including the effect of temperature and pressure of 
the reactive environment. In this letter, we address the 
issue of stability and metastability using our model sys¬ 
tem that has significance in many practical applications. 
Our model includes free metal (Li) clusters in a (a) ni¬ 
trogen only atmosphere and then in (b) an atmosphere 
of both nitrogen and hydrogen. It should be noted that 
we have addressed the situation under thermodynamic 
equilibrium, while the real system might not be always 
in thermodynamic equilibrium. But we emphasize that, 
a thermodynamic phase diagram acts as guideline and 
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FIG. 1: Phase diagram for LisNj; clusters in an nitrogen at¬ 
mosphere. The geometries are optimized with PBE-|-vdW 
and the total energies are calculated with PBEO-l-vdW. 

depicts its limit for predicting properties and functions 
of real materials. 

We have considered first a wide range of Li 3 ;Nj^ clusters, 
where a;=l, 2,..., 6, 9, 12, 15, 18 and y is determined by 
thermal equilibrium with the environment at given tem¬ 
perature T and partial nitrogen pressure pn 2 [IZ] ■ To get 
the minimum energy conhgurations, for each stoichiom¬ 
etry, the total energy is minimized with respect to both 
geometry and spin state. Very unexpectedly, our results 
reveal that though in bulk LiaN is the most stable phase, 
in small clusters, the stoichiometry of the most stable 
phase is quite different. 

The low energy structures (including the global mini¬ 
mum (GM)) are generated from an exhaustive scanning 
of the potential energy surface (PES) from our recent 
implementation of cascade genetic algorithm. The term 
cascade means a multi-stepped algorithm where succes¬ 
sive steps employ higher level of theory and each of the 
next level takes information obtained at its immediate 
lower level. Typically the cascade GA starts from a crude 
but sophisticated classical force field and goes up to den¬ 
sity functional theory using hybrid functionals. This GA 
algorithm’s implementation is thoroughly benchmarked 
and it’s efficiency is validated (w.r.t more advanced the¬ 
ory) in details in Ref [571ISH] ■ 

We have performed the density functional theory (DFT) 
calculations using FHI-aims, which is an all electron 
code with numerical atom centred basis sets [39]. The 
low energy GA structures are further optimized at a 
higher level settings, where energy minimization is per¬ 
formed with vdW-corrected-PBE-|-vdW SO] functional, 
“tight - tier 2” settings, and force tolerance was set 
to better than 10“^ eV/A. The van der Walls correc¬ 
tion is calculated as implemented in Tkatchenko-Scheffier 
scheme m- Finally the total single point energy is 
calculated afterwards on top of this optimized struc¬ 


ture via vdW-corrected-PBEO [15] hybrid xc functional 
(PBEO-l-vdW), with “tight - tier 2” settings. [48| We 
find that PBE-|-vdW strongly overestimates stability of 
clusters with larger y values in Li 2 ,Ny clusters. This re¬ 
sults a qualitatively wrong prediction that adsorption of 
N 2 could be favored over desorption up to a large ex¬ 
cess of nitrogen. Such behavior is not confirmed by hy¬ 
brid functional [e.g HSE06 [43], PBEO] as employed in 
our calculations. The difference in energetics of PBEO 
and HSE06 is always within 0.05 eV. The spin states 
of the clusters are also different as found by PBE and 
PBE0/HSE06. In view of this, all our results are thor¬ 
oughly tested and benchmarked w.r.t hybrid functionals 
(PBEO) using “tight” numerical settings and tier 2 basis 
set [55] . 

The free energy of the low energy isomers [55] (in the 
PES) is then calculated as a function of T and pns for 
each stoichiometry using the ab initio atomistic thermo¬ 
dynamics (aiAT) approach. The concept of aiAT is ear¬ 
lier developed and successfully applied initially for bulk 
semiconductors [H], and later applied to the study of 
oxide formation at the surface of some transition metals 
and other materials [JS] . We have recently extended this 
approach to cluster systems m following our detailed de¬ 
scription as in Ref [38]. Therefore, from different cluster 
compositions and structures with the lowest free energy 
the thermodynamic phase diagram can be constructed 
as a function of T and pn 2 ■ One such phase diagram is 
shown in Fig. At a given T, pn 2 i and y, the stable 
stoichiometry of a Li^jN^ cluster is determined via aiAT, 
i.e., by minimizing the Gibbs’ free energy of formation 
^0'/(r, PNa)- 

^G'/(r,pN2) = Fu^NyiT)-Fi,i^{T)-yxy]si{T,p^^) (1) 

Here, TLi„,N (T) and Fi,i^ (T) are the Helmholtz free en¬ 
ergies of the LiajNj, and the pristine Lia, cluster [50] and 
/iN(T,pN 2 ) is the chemical potential of nitrogen. As ex¬ 
plained in Ref [55], (T) and (T) are calculated 

using DFT information and are calculated from the sum 
of DFT total energy, DFT vibrational free energy in the 
harmonic approximation, as well as translational, rota¬ 
tional, symmetry and spin-degeneracy free-energy con¬ 
tributions. The dependence of on T and pjq^ 

is calculated using the ideal (diatomic) gas approxima¬ 
tion with the same DFT functional as for the clusters. 
The phase diagram for a particular Li 2 ;Ny is constructed 
by identifying the lowest free-energy structures at each 
T, pn 2 - A.S a representative example, we show in Fig. [^ 
the phase diagram for x = 3. Note that surprisingly in 
small cluster, we find that at a realistic T and pn 2 > N-rich 
phases are becoming more stable over the conventionally 
known most stable phase LiaN as in bulk. We have then 
extended our observation for a wide range of clusters of 
Li 2 ,Ny conhguration varying x as 1, 2,.., 6, 9, 12, 15, 18 
and y is increased as 1, 2, 3,.. as per the highest possible 
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FIG. 2: (color online) The most stable Lia^N^ clusters at var¬ 
ious temperatures and pressures under thermodynamic equi¬ 
librium. The geometries are optimized with PBE-|-vdW, and 
the electronic energy is calculated using PBE0-|-vdW. 


value at each x in the thermodynamic condition as em¬ 
ployed in the phase diagrams. For example, for a: = 3, we 
see from Fig.[2the highest possible y is 11. We find that 
the N-rich phase in these small clusters are consistently 
becoming quite stable over a wide range of en¬ 

vironmental conditions. This is shown in Fig. where 
we have plotted the most stable stoichiometry of y at 
different with varying x values. The black line 

represents the usual Li:N=3:l stoichiometry as in bulk. 
Note that most of the stable phases are above this black 
line and are significantly high N-rich especially for small 
clusters. As expected, due to quantum size effect, the 
reactivity of the larger clusters get decreased with grow¬ 
ing size. But still even at larger sizes, the clusters are 
more N-rich as compared to their bulk limit. Only at 
high temperature (> 600 K) and low pressure (< 10“® 
atm) the clusters are N-deficient. 

Subsequently, we explored the efficiency of the N-rich 
phases of these cluster as hydrogen storage material. 
However, the procedure for this is not easy as one has 
to understand under which temperature conditions and 
in which N 2 and H 2 (partial) pressure range, these nan¬ 
oclusters become more stable, i.e. under which optimum 
conditions the N 2 and H 2 adsorption can take place on 
the surface of a given cluster size. Note that under re¬ 
action conditions, the nano-cluster comprises of a wide 
range of structures including different number of atoms 
with various oxidation states, all of which could be ac¬ 
tive to some extent in the reaction. Some obvious ques¬ 
tions arise naturally, for example, “which are the species 
present in the real reaction and what are their struc¬ 
tures?”. It’s also interesting to understand “how the 
nano-cluster changes their structure and properties upon 
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FIG. 3: (color online) Phase diagram of LiN^Hz at varying 
pressure of ph 2 and pn 2 of two different temperatures T = 
300 K (top) and T = 600 K (bottom) 


adsorption of a different ligand molecules? We have ad¬ 
dressed these question on considering another ligand (e.g 
H 2 ) adsorbed onto LiNy clusters forming LiN^H^ clusters 
and tried to understand the stability of the entire system 
in presence of both N 2 +H 2 pressures. In order to ful¬ 
fil these aims, small clusters LiNj^H^ are first thoroughly 
scanned from our cascade genetic algorithm implementa¬ 
tion and their stability is studied using atomistic thermo¬ 
dynamics to mimic an atmosphere composed of N 2 and 
H 2 gases at realistic temperature and pressure conditions. 
As depicted in Eq.[^ here to address the (meta-)stability 
of LiNyH^ clusters, in presence of one more additional 
ligands i.e. H 2 , the equation gets modified as below: 


^G'/(T,pn 2,TH2) = - ALi(T) 

-y X /xn(T,pn 2 ) - -2 X Mh(7’,PH2) 

Here, (T’) and Fi^i(T) are the Helmholtz free 

energies of the LiN^Hj, and the pristine Li cluster (in 
this case single Li atom) respectively. Like before as in 
Eq.[TJ here also we have calculated this from the sum 
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FIG. 4: (color online) Release of three units of ammonia from 
LiN 4 Hii before the final product LiNH 2 is formed at a varying 
T , PNa , PHa • 


of PBEO DFT total energy, DFT vibrational free energy 
(in the harmonic approximation), and translational, ro¬ 
tational, symmetry, spin-degeneracy free-energy contri¬ 
butions. /X 7 v(T,PN a) are the chemical po¬ 

tential of nitrogen and hydrogen respectively. 

In Fig. we have compared the phase diagrams of the 
stability of LiN^Hz clusters, varying H 2 and N 2 pressures, 
at two different temperatures (viz. 300 K (top) and 600 
K (bottom)). It’s clear from Fig. |^that at T=300 K, the 
small clusters can adsorb quite a substantial amount of 
hydrogen at a realistic temperature and pressure. Note 
that in bulk LisN, the highest amount of hydrogen up¬ 
take is possible when it forms lithium amides/imides 
(LiNH 2 /Li 2 NH). But in small clusters, the situation is 
very different as LiN 4 Hii is thermodynamically the most 
stable cluster and the material is stable over a wide ex¬ 
perimentally achievable region in the phase diagram at 
T=300 K. It should be noted that for an efficient hy¬ 
drogen storage materials, the desorption kinetics should 
also needs to be optimal. We have provided a similar 
phases diagram in the bottom panel of Fig. at T = 
600 K. We find that at this temperature first the stabil¬ 
ity region of LiN 4 Hii is reduced, and less H-rich phases 
(viz. LiNsHg, LiN 2 H 5 , LiNH 2 ) become more stable. This 
means the cluster tends to leave hydrogen at higher tem¬ 
perature, while the H-rich phase being absolutely stable 
at room temperature. On a critical look to the phase dia¬ 
gram, it appears to be that under a fixed pressure region 
of both PH 2 and PN 2 ) if simply increase temperature 
these type of clusters have a tendency to release ammo¬ 
nia (NH 3 ) molecule. This is shown more clearly in Fig.|^ 
It’s evident from Fig. |^that on increasing temperature 


these nano-clusters will form a reversible cycle under a 
constant ph 2 and PN 21 where the following reaction is 
going on as shown in Fig. 

LiN4Hii ^ LiNaHg-fNHg ^ LiN2H5-tNH3 ^ LiNH2-fNH3 

( 3 ) 

Clearly, just by changing temperature we can get three 
units of NH 3 out of this LiN 4 Hii. Therefore, getting such 
an automatic release of ammonia solves the problem of 
hydrogen storage, as hydrogen can be released econom¬ 
ically from ammonia on-demand, without the need for 
high-pressure or cryogenic storage |46) . 

In summary, we have employed our massively parallel 
cascade genetic algorithm to scan the potential energy 
surface of a large number of Li^^Ny clusters. We find 
that a local/semi-local functional is not correct even to 
find a reasonably qualitative prediction. In view of this, 
we have employed more advanced hybrid DFT functional 
PBFO to calculate accurately the energetics throughout 
our calculations. We find by applying aiAT, that the 
behaviour of small Lij^Ny clusters are very different than 
their bulk, which may be due to quantum confinement ef¬ 
fect. In small clusters, over a wide range of sizes, N-rich 
phases are consistently stable at thermodynamic equilib¬ 
rium under experimentally achievable temperature and 
pressure pn 2 - We, therefore, presumably for the first 
time introduce a new class of clusters, which has been 
overlooked in the past. We went one step further to un¬ 
derstand the effectiveness of these clusters in the field 
of hydrogen storage. We generated a large number of 
Li 2 ,NyH 2 clusters out of which one of the test cases is 
shown here for LiN^H^. These clusters not only can store 
significantly high amount of hydrogen but also possesses 
a fast release kinetics of hydrogen. Therefore, these ma¬ 
terials should have pretty wide application in the field of 
hydrogen storage. 

*Please address all correspondence to: Saswata Bhat- 
tacharya (saswata@physics.iitd.ac.in) 
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